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Abstract: We investigate the relativistic equation of state of hadronic matter and quark-gluon plasma at finite 

temperature and baryon density in the framework of the non-extensive statistical mechanics, characterized 
by power-law quantum distributions. Wc impose the Gibbs conditions on the global conservation of baryon 
number, electric charge and strangeness number. For the hadronic phase, we study an extended relativistic 
mean-field theoretical model with the inclusion of strange particles (hyperons and mesons). For the quark 
sector, we employ an extended MIT-Bag model. In this context wc focus on the relevance of non-extensive 
effects in the presence of strange matter. 
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1. Introduction 

The behavior of hot, dense nuclear matter is usually described in terms of the so-called relativistic mean field 
models, wherein thermodynamical quantities are usually obtained through the common Boltzmann-Gibbs statis- 
tics. But, in recent years, there have been growing theoretical and experimental indications that, under extreme 
conditions and during the possible phase transition to quark-gluon matter reachable in high energy heavy ion col- 
lisions, strong dynamical correlations, long-range color interactions and microscopic memory effects can take place 
[1-18]. In this context the non-extensive statistical mechanics proposed by Tsallis [1-3] can be used to describe 
and investigate such physical phenomena. Nonextensive statistical effects should strongly affect the finite tem- 
perature and nuclear density Equation of State (EOS) [19-23]. In fact, by varying temperature and density, the 
EOS reflects (in terms of the macroscopic thermodynamical variables) the microscopic interactions between the 
different nuclear matter phases. This means that the physical proprieties of the system can be sensibly changed 
with respect to the standard Boltzmann-Gibbs statistics. 
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Unfortunately, extracting information about the EOS at different densities and temperatures by means of heavy 
ion collisions is very difficult and can be realized only indirectly by comparing the experimental data with different 
theoretical models such as, for example, fluid-dynamical models [24]. ft is also relevant that a non-extensive 
hydrodynamic model for multiparticle production processes has been proposed [25]. 

In this paper, we extend previous investigation (Ref. [22]) to the regime of heavy- ion collisions, by by considering 
the physical conditions in which strange (hadronic and quark) particles can take place, and we analyze the mixed 
phase region following the Gibbs condition for the phase equilibrium by requiring the global conservation of the 
three charges: baryon number, electric charge and zero net strangeness. 

We show that, in the presence of non-extensive statistical effects, there is a remarkable variation in many physical 
features such as a strong enhancement of the strange hadronic particles production both in the hadronic and in 
the mixed phase, together with a significant variation in the particle concentrations and a softening of the pressure 
in the mixed phase. 

2. Nonextensive hadronic equation of state 

Nonextensive statistical mechanics, introduced by Tsallis [1-3], is a generalization of the common Boltzmann- 
Gibbs statistics, based upon the introduction of the following entropy 

sM = ^y (i- J If Wan) , (| /(x) dft = 1) , (1) 

where /(x) is a normalized probability distribution, with x and dQ, denoting, respectively, a generic point and the 
volume element in the corresponding phase space. The generalized entropy has the usual properties of positivity, 
equiprobability, concavity and irreversibility, and in the limit of q — > 1, the entropic form (1) becomes additive 
and reduces to the standard Boltzmann-Gibbs entropy 



Si = - [f(x) ln/(x)<ifi. (2) 

A second crucial assumption on non-extensive statistics is the introduction of the q-mean value (or escort mean 
value) of a physical observable j4(x) 

/A(x) [/(x)]«dn 

{A} «- nn^dn ■ (3) 

The non-extensive statistical effects vanish approaching zero temperature, and such a formalism can be considered 
as an appropriate basis to deal with physical phenomena in which strong dynamical correlations, long-range 
interactions and microscopic memory effects take place [1-4]. 

The probability distribution can be obtained maximizing the measure S q under appropriate constraints related 
to the previous definition of the q-mean value. In this context, it is important to observe that the Tsallis classical 
distribution can be seen as a superposition of Boltzmann distributions with different temperatures whose mean 
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values correspond to the temperature in the Tsallis distribution. The non-extensive q parameter is related to 
the temperature fluctuation and describes the spread around the average value of the Boltzmann temperature 
[4]. Moreover, let us remember that, in the diffusional approximation, a value q 7^ 1 implies anomalous diffusion 
among the constituent particles with superdiffusion if q > 1, and subdiffusion if q < 1 [26]. 

From the above, we can obtain the associate quantum mean occupation number of particles species i in a grand 
canonical ensemble. For a dilute gas of particles and for small deviations from the standard statistics (9 » 1), the 
occupation number can be written as [19, 20] 

(4) 



e q (P(Ei - m)) ± 1 



where /3 = 1/T and the sign option (±1) is for fermions and bosons respectively. Furthermore, in Eq.(4), following 
Ref. [20], for q > 1, we have 

[l + (g-l)a:] 1/(9-1) ifx>0; 



e q (x) 
whereas, for q < 1, 



e q (x) = < 



(5) 

[l + (l-q)x] 1/{1 - q) ifx<0, 



[1 + (q- l)x] 1/{q - 1) ifx<0; 

(6) 

[l + (l-q)x] 1/{1 - q) ]£x>0. 



Naturally, as q — > 1, the above quantum distribution reduces to the standard Fermi-Dirac and Bose-Einstein 
distributions. 

Let us observe that when the entropic q parameter is smaller than one, the high-energy tail of the above particle 
distribution is depleted; whereas when q is greater than one, it is enhanced. Hence the non-extensive statistics 
entails a sensible difference of the power-law particle distribution shape in the high energy region compared with 
the standard statistics. From a phenomenological point of view, the non-extensive index q is considered here 
as a free parameter, even if it is actually not because it should depend, in principle, on the physical conditions 
generated in the reaction and on the fluctuation of the temperature, and also be related to microscopic quantities 
(such as the mean inter-particle interaction length, the screening length and the collision frequency into the parton 
plasma) . 

In the present investigation we are going to study small deviations from the standard statistics and values of q > 1 
only, because these values were obtained in several phenomenological studies of high energy heavy ion collisions 
(see, for example, Ref.s [12, 15, 18, 27]). 

In this framework, we analyze hadron interaction through the relativistic mean field (RMF) model [28-31], in 
order to investigate the equation of state of nuclear matter at finite baryon density and at temperatures reachable 
in high energy heavy ion collisions. 
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The Lagrangian density for the full octet of baryons (p, n, A, E + , E°, E , H°, S ) can be written as 

£octet = ^2 Tpi [i 7 M <9 m — (Mi — g a i a) — g^i 7 M ^ — g P i 7 M t- p M ] tpi + - (d^ad^a — mla 2 ) — U (a) 

i 

+ \ m l> <*V wM + i c (gl N oj m cj m ) 2 + i m 2 p fa ■ p M - jF^F^ - ^G^&" , (7) 

where the sum runs over all baryon octets, Mi is the vacuum baryon mass of index i, t denotes the isospin operator 
which acts on the baryon and U(a) is the nonlinear self-interaction potential of a meson 

U(o) = ^a(gaNcr) 3 + i fe (,g CT jv ct 4 ) , (8) 

introduced by Boguta and Bodmer [29] in order to achieve a reasonable compressibility at the saturated density 
of nuclear matter. 

The meson fields are coupled with the baryon octet through opportune model-dependent coupling constants. In 
this scheme, the effective baryon mass is given by M* = Mi — geno, and the scalar and vector baryon density pf 
and pf are given, respectively, by 



(9) 



P?=H J |4[»ii(l!)-ni(«;)], (10) 

where ji = 2Ji + 1 is the degeneracy spin factor and ru(k) and ru(fc) are the g-deformed particle and antiparticle 
distributions function given in Eq.s(4)-(6); for example for q > 1 and /3(E* — \p*\) > 0, we have 

Mk) = [i + (g-i)^r(fe)+M:)]V( 9 -i) + i ■ ^ 



The baryon effective energy is defined as E* (k) = \/k 2 + Mi* 2 and the effective chemical potentials p* are given 
in terms of the meson fields as follows 

Pi = Pi - g^jiu - g P %T-ap , (13) 

where pi are the thermodynamical chemical potentials pi — de/dpi. At zero temperature they reduce to the 
Fermi energies and the non-extensive statistical effects disappear. The meson fields are obtained as a solution of 
the field equations in mean field approximation and the related meson-nucleon couplings constant (g a N, gu>N and 
<7 p jv) will be fixed to the parameters set marked as GM3 of Ref.[30]. 

Because we are going to describe a finite temperature and density nuclear matter with respect to strong interaction, 
we have to require the conservation of three 'charges': baryon number (B), electric charge (C) and strangeness 
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number (5). For this reason the system is described by three independent chemical potentials: hb, \ac and fj,s, 
respectively, the baryon, the electric charge and the strangeness chemical potentials. Therefore, the chemical 
potential of particle of index i can be written as 

Mi =bi(iB + a (J-c + Sifis , (14) 

where bi, a and Sj are, respectively, the baryon number, the electric charge and the strangeness quantum numbers 
of the i-th hadronic species (baryons and mesons). 

The thermodynamical quantities can be obtained from the baryon grand potential Qb in the standard way. More 
explicitly, the baryon pressure Pb = —Qb/V and the energy density can be written as 

Pb = ^J]^ J -j^y -^j^ +n q i(k)] - ^mla 2 - U(a) + ^mtu; 2 

+ ^c(g^ N bj) 4 + ^m 2 p p 2 , (15) 

en = X>/ ^m(k)[nUk) + n«(k)] + ±mla 2 + U(*) + ±rniu 2 

+ ^c(g^NUj) 4 + ^m 2 p p 2 . (16) 

Moreover, as temperature increases and baryon density decreases, the mesonic contribution to the total thermo- 
dynamical potential becomes increasingly important. For simplicity, we treat the mesons as an ideal Bose gas so 
we evaluate the pressure Pm , energy density €m , particle density p M of mesons as 

^ = sE-»/|$r^W) + W], (17) 



3 

M 



Iwf ^(fc)K(fc) + fi«(fc)], (is) 



d 3 k 



Pi =1iJ 7^)3 M*) -**(*)]. (19) 

where 7j = 2Jj + 1 is the degeneracy spin factor of the j'-th meson and the functions rij(k) fij(k) are the g- 
deformed boson particle (antiparticle) distributions of the j'-th meson. The total pressure and energy density is 
given as usual by P tot = P B + Pm and e to t = e B + €m- 

3. Quark-gluon equation of state 

In this work we use a simple effective MIT bag model to describe the quark phase. All the non-perturbative effects 
are simulated by the bag constant B which represents the pressure of the vacuum. It is well known that, using 
the simplest version of the MIT bag model, at moderate temperatures the deconfinement transition takes place at 
very large densities if the bag pressure B is fixed to reproduce the critical temperature computed in lattice QCD. 
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On the other hand, there are strong theoretical indications that at moderate and large densities (and not too 
large temperatures) diquark condensates can form, whose effect can be approximately accounted for by reducing 
the effective bag constant [32]. As proposed in Ref. [33], a phenomenological approach can therefore be based on 
an effective "bag constant" B e g depending on the baryon chemical potential, ft can be written as 

Bcff = (Bo - 5oc)/(l + exp[( A t s - po)/o]) + Boo , (20) 

where we have set B^ 4 = 250 MeV (bag constant at vanishing ^b), BiJ, 4 — 160 MeV (bag constant at very large 
/is), /Ho = 600 MeV and a = 320 MeV. The above parameter values have been fixed by requiring that, at low hb, 
the critical temperature is « f70 MeV for q = 1, while the other constraint is the requirement that the mixed 
phase starts forming at a density slightly exceeding 3po for a temperature of the order of T m 90 -j- 100 MeV (as 
also suggested, e.g., by [34]). 

Following this line, the pressure, energy density and baryon number density for a relativistic Fermi gas of quarks 
can be written, respectively, as 

^=yE/ o ^s^H(k)+rf f (k)]-B eS , (21) 
e = IS E J o §^3 e f K (A) + n} (*)] + B cS , (22) 

p-yE( ^ [«/(*) -»/(*)]. (23) 

where the quark degeneracy for each flavor (/ = u,d,s) is 7/ = 6, rif(k) and nj{k) are the g-deformed particle 
antiparticle quark distributions. Gluons and light quarks (u, d) are considered as massless point-like particles, 
while for strange quarks (s,s) we consider a finite mass m s = f50 MeV. 

4. Mixed hadron-quark-gluon phase 

In this section, we investigate the hadron-quark phase transition at finite temperature and baryon chemical 
potential. In a theory with only giuons and no quarks, the transition turns out to be of first order. Since the u 
and d quarks have a small bare mass, while the strange quark has a somewhat larger mass, the phase transition 
is predicted to be a smooth cross-over. In fact, various results from QCD-inspired models indicate that, as 
the baryon chemical potential increases in the phase diagram, a region of non-singular but rapid crossover of 
thermodynamic observable around a quasi-critical temperature, leads to a critical endpoint (CEP), beyond which 
the system shows a first-order phase transition from confined to deconfined matter. However, the existence or 
exclusion of a CEP has not yet been well confirmed by QCD lattice simulations. Since it occurs over a very 
narrow range of temperatures, the transition, for several practical purposes, can still be considered of first order. 
Indeed the lattice data with 2 or 3 dynamical flavours are not precise enough to unambiguously disentangle the 
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difference between the two situations. Thus, by considering the deconfinement transition at finite density as a 
first-order one, a mixed phase can be formed, which is typicaily described using the two separate equations of 
state: one for the hadronic phase, and one for the quark phase. 

To describe the mixed phase, we apply the Gibbs formalism to systems with more than one conserved charge 
[35, 36], requiring that baryon number, electric charge and strangeness number are preserved. The main result is 
that, at variance with the so-called Maxwell construction, the pressure in the mixed phase is not strictly constant 
and therefore, for instance, the nuclear incompressibility does not vanish. 

The structure of the mixed phase is obtained by imposing the Gibbs conditions for chemical potentials (/x^' = 
\ Mc^ = Pc \ Ps = Ps^) an d pressure (P H — P®) and by requiring global conservation of the total 
baryon, electric charge and strangeness densities in the hadronic (H) and quark (Q) phases: 

P (H) (T, n B , fi c , Ps) = P (Q) (T, p B ,pc,Ps) , 

PB = 0-~ X)pB + XPi , 

PC = (1 - x)pc +XPc > 

PS = (1 - x)Ps + XP Q s , (24) 

where \ i s the fraction of quark-gluon matter in the mixed phase. In this way we can find the phase coexistence 
region in the (T, pb , pc , ps) space. At fixed T and [is, the charge pc and strangeness ps chemical potentials 
are obtained by fixing the total electric charge Z/A (for example, Z/A = 0.4 for lead-lead heavy ion collisions) 
and the total strangeness neutrality by the conditions 

pc = -jPb, Ps — ■ (25) 

In Fig. 1, we report the variation of the pressure as a function of baryon density (in units of nuclear saturation 
density po — 0.153 fm~ 3 ) by considering the degrees of freedom of strange matter in the hadronic and quark 
phase, for different values of q. The temperature is fixed at T — 120 MeV and the electric charge fraction to 
Z/A — 0.4. In presence of non-extensive effects the pressure results to be considerably increased even for small 
deviations from the extensive standard statistics. 

It is interesting to observe that pressure (as a function of baryon density) is stiffer in the pure hadronic phase, but 
a strong softening appears in the mixed phase. This feature is remarkable evident by increasing the non-extensive 
entropic parameter q; it implies an abrupt variation in the incompressibility and may be particularly important in 
identifying the presence of non-extensive effects in high-energy compressed baryonic matter experiments. Indirect 
indications of a sensible softening of the EOS at the energies reached at AGS have already been raised elsewhere 
[24, 33, 37-39]. 

In Fig. 2, we show the phase diagram in the pB — T plane for different values of q. The curves labeled with p\ r and 
pir denote, respectively, the beginning and the end of the mixed phase. In presence of non-extensive statistical 
effects, we observe a remarkable lowering of the critical maximum temperature at vanishing baryon density ps- 
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Figure 1. Pressure as a function of the baryon density for different values of non-extensive parameter q (GM3 parameters set). 
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Figure 2. Phase diagram in the ps — T plane for different values of q. The curves labeled with p l CT and p^r stand for the beginning 
and the end of the mixed phase, respectively. 

In Fig. 3, we show the most relevant particle concentration Yi as a function of the baryon density in the pure 
hadron phase, mixed phase and quark phase, at fixed temperature T — 120 MeV for q — 1 (left panel) and 
q = 1.1 (right panel). As we can see, in presence of non-extensive statistical effects, the particle concentrations 
clearly differ. In particular, we observe a strong reduction of the neutron and proton fractions with a consistent 
enhancement of the hyperon and meson fractions, whereas the quark concentrations are not sensibly modified. 
Finally, in Fig. 4, we show the net strangeness density p s for the different particle species (baryons, mesons and 
strange quarks) as a function of baryon density at fixed temperature T — 120 MeV and for different values of q. 
In presence of non-extensive statistical effects (q — 1.1), we observe a strong enhancement of the strange hadronic 
particles production, especially in the pure hadronic phase. At the beginning of the mixed phase (ps = 2po), we 
can observe a maximum in the baryon and mesons net strangeness density, of about 4 times the values obtained 
in the standard statistical case (q — 1), after which it decreases until zero at the beginning of the quark phase. 
Nonextensive statistical effects have instead a lower impact on the quark sector and the variation of the quark 
strange density results to be negligible. 
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Figure 3. Particle concentrations as a function of the baryon density in the pure hadronic phase, mixed phase and quark-gluon 
matter for q = 1 (left panel) and q = 1.1 (left panel), T = 120 MeV and Z/A = 0.4. 
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Figure 4. Net strangeness densities of baryons (B), mesons (M) and strange quarks (s) for different values of the non-extensive 
parameter q at a fixed temperature of T — 120 MeV, in pure hadronic phase, mixed phase and quark-gluon phase. 



To better understand the enhancement of the strange hadronic particle in presence of non-extensive statistical 
effects observed in the previous two figures, in Fig. 5 we report the effective mass M* (with respect to the 
corresponding vacuum mass Mi) for the most relevant baryons (nucleons N, A, E) as a function of the baryon 
density at the same values of q and temperature reported in Fig.s 3 and 4. 

The behavior of the effective baryon mass as a function of the baryon density ps (or baryon chemical potential 
Pb) is related to the behavior of the a field by means of the relation M* = Mi — g a tcr (see Section 2). According to 
the results of Ref. [22] , we find that at low pB , in presence of non-extensive effects, the value of the a meson field 
is increased respect to the standard case and decreased at higher values of ps- This important feature is due to 
the fact that for q > 1 and fixed baryon density (or /is), the (normalized) mean occupation function is enhanced 
at high values of its argument and depressed at low values. Being the argument of the mean occupation function 
Xi — P{E* — p*), in the integration over momentum (energy), at lower fiB (corresponding to lower values of the 
effective particle chemical potential p,*) the enhanced Tsallis high energy tail weighs much more that at higher ps 
where depressed low energy effects prevail and the mean occupation number results to be bigger for the standard 
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Figure 5. Effective masses ratios (M* /Mi) of the most relevant baryon particles (nucleons N, A, S) versus the baryon density 
Pb, for different values of the non-extensive parameter q at a fixed temperature T — 120 McV. 

Fermi-Dirac statistics. Concerning the antiparticle contribution, the argument of Hi is Xi = P{E* + p*) and the 
Tsallis enhancement at high energy tail is favored also at higher p,s- At the same time, higher temperatures 
(where antiparticle contribution are more relevant) reduce the value of the argument of rii and Hi, favoring the 
extensive distribution. These effects are much more evident for the scalar density ps (self-consistently related to 
the a meson field) where appears (rii) q and particle and antiparticle contributions are summed. As a consequence, 
the baryon effective masses become, respect to the standard case, smaller for lower values of ps and bigger for 
higher values. As can be seen in Fig. 5, this feature appears to be more significant as a percentage for lighter 
baryons, therefore, at large baryon densities {ps ~ 1 -5- 2po), the nucleons effective mass is enhanced respect to 
the standard statistical case (q = 1) as a sensible greater percentage with respect to the A and £ effective masses. 
This matter of fact favors the formation of hyperons respect to nucleons at a fixed baryon density and in presence 
of non-extensive statistical effects. Moreover, being the strangeness number globally conserved, an increase of 
hyperon particles implies, in the pure hadron phase, a corresponding increase of strange meson particles in order 
to satisfy the condition of zero net strangeness. This effect is also present in the mixed hadron-quark-gluon phase 
because of the strange quark formation is not sensible affected by non-extensive statistical effects (such a behavior 
could be due to the simple MIT bag model adopted in this paper; in this sense it could be very interesting future 
investigations with more sophisticated SU(3) chiral quark models). 

In light of the above results, a significative enhancement of the strange to non-strange particle ratios at finite 
temperature and baryon density could be an evident signature of the relevance of non-extensive statistical effects 
in the future compressed baryonic matter experiments. 
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